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Abstract 

Hadronic correlation functions at finite temperature in QCD, with four flavours 
of dynamical quarks, have been analyzed both above and below the chiral symme- 
try restoration temperature. We have used both point and extended sources for 
spatial as well as temporal correlators. The effect of periodic temporal boundary 
conditions for the valence quarks on the spatial meson correlators has also been in- 
vestigated. All our results are consistent with the existence of individual quarks at 
high temperatures. A measurement of the residual interaction between the quarks 
is presented. 



1 Introduction 



The spatial correlation functions of operators with hadronic quantum numbers yield 
screening masses for static excitations in the QCD plasma ||. These give informa- 
tion on the physical excitations and interactions in QCD at finite temperatures. 
Chiral symmetry restoration at temperatures T > T c is signalled by the degeneracy 
of screening lengths obtained from pairs of opposite parity channels [|l], 0, [3|. At 
such high temperatures, the vector (V), as well as the pseudovector (PV), screening 
mass is nearly twice the lowest Matsubara frequency, Q = ttT, whereas the baryon 
screening mass is close to thrice this value. This indicates that correlations in these 
channels are mediated, respectively, by the exchange of two and three weakly inter- 
acting quarks ||. 

Close to T c , the scalar (S) and the pseudoscalar (PS) 'meson' screening mass is 
significantly smaller than the V/PV screening mass, but approaches the latter when 
the temperature is raised to 2T C [|[. Thus, at high temperature, this channel is also 
correlated through the exchange of two weakly interacting quarks. Nevertheless, 
close to T c , such an interpretation does not seem to hold. 

Since the lowest mass hadron is unlikely to change from a mesonic collective state 
to a quark-like quasi-particle at a non-critical point, a pleasantly consistent picture 
of the excitation spectrum would be obtained if, above T c , the S/PS correlations 
could be shown to be due to the propagation of two unbound quarks, possibly still 
strongly interacting close to T c . A means of checking this was suggested and applied 
to quenched QCD [ffl . This method also yielded measurements of effective couplings 
in different spin channels. It was seen that such a coupling is indeed larger in the 
S/PS channel than in the V/PV channel. In this paper we extend such studies to 
QCD with four flavours of dynamical staggered quarks. 

Hadronic correlators, and similarly the spatial Wilson loops || and spatial four 
point functions with point-split sources ||, should provide information on the non- 
perturbative long distance structure of the quark gluon plasma. While the spatial 
and temporal quark propagators yield identical information about quark screening 
masses |7j], this seems to be in conflict with recent studies of temporal and spatial 
hadron propagators M. In the latter case larger screening masses were extracted 
from temporal than from spatial correlators, although one would have expected 
to find smaller values as the lowest Matsubara mode is zero in this case. The 
determination of low momentum excitations from temporal propagators at finite 
temperature is, however, difficult, due to the short "time" extent of the lattice 
(0 < t < 1/T). This leads to a superposition of many high momentum excitations. 
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We will try to clarify this situation here by projecting onto low momentum modes, 
by using wall source operators. 

Furthermore, we will study the sensitivity of the extracted spatial screening 
lengths to changes of the temporal quark boundary conditions. While bosonic modes 
(bound states in mesonic operators) will not be sensitive to changes of the usual anti- 
periodic boundary conditions to periodic ones, fermionic modes (freely propagating 
quarks) will be sensitive to these changes. We will also extend the studies of scalar 
and vector channel couplings |4[] to the case of QCD with dynamical fermions. 

The conclusion from these investigations is that temporal and spatial propagators 
yield identical information on the spectrum of finite temperature QCD below and 
above the chiral phase transition. The behaviour of the (pseudo-)scalar screening 
length suggests that there are no bound states in this quantum number channel 
above T c . 

This paper is organized as follows. In section 2 we introduce the observables 
studied, and describe their behaviour for free staggered fermions, i.e., in the lowest 
order of perturbation theory. Details of our measurements and the results are given 
in section 3. A summary and conclusions are presented in the section 4. 

2 Spatial and Thermal Direction Correlators 

In this section we consider correlation functions between operators separated either 
in the spatial or the thermal direction. They are constructed on lattices of size 
V = N T x N% with N a > N T . Unless otherwise mentioned, we shall use the lattice 
spacing as the unit of length. The shorter direction introduces the temperature, via 
N T = 1/T. Spatial and thermal direction correlation functions with given external 
meson momentum P M are defined as 

Gf(x 3 ,P) = -^—J2(H(x 3 ,x)H^0,6))e^(iP-x), 
a T x 

G*(x ,P) = ^£(P(x , *)#+((), 0)>exp(*P.*). (2.1) 

For any 4-momentum P M , we shall use the notation P to denote the 3-momentum 
(Po, Pi, P2), and P for (Pi, P 2 , P3). Here H denotes an operator carrying mesonic 
quantum numbers. By an abuse of language, we shall call G s and Gt 'meson cor- 
relators'. The name is not supposed to carry any prejudice about the existence or 
absence of mesons in the theory. 

In perturbation theory, these correlation functions can be described in terms of 
the exchange of two (interacting) fermions. We set out our notation and introduce 
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some features typical of non-interacting staggered fermions below. This corresponds 
to the leading order, O(g ), of perturbation theory. In section 3 the behaviour of 
the measured correlation functions will be compared with these features. 

2.1 Spectral representations 

Hadronic correlation functions at T = are used to extract the mass of the lowest 
excitation with given quantum numbers. A spectral representation of correlation 
functions with quantum numbers H, in the form 

G H (t) = £Ae- m K (2.2) 

i 

can be used to show that this necessitates the measurement of G H (t) at Euclidean 
time separations much larger than the splitting between the lowest and the first 
excited states, i.e., t 3> 1/ — m Q H ). In the absence of any a priori knowledge of 
this splitting, it is necessary to take as large a time separation as possible, and to 
perform checks on the estimate of m? H so obtained. All this is well known. 

At finite temperatures the physical extent of the thermal direction, r, is bounded 
in physical units to r < 1/T. Due to the periodicity (or antiperiodicity) imposed 
to obtain the thermal ensemble, correlation functions can be followed only to r < 
1/(2T). In the general case this may not allow the extraction of the lowest lying state 
in the manner discussed above. One then has to use either a complicated ansatz for 
the correlation function, involving several exponentials to approximate eq. ( j2.2| ) ||, 
or a more complicated operator which projects out only the low momentum modes. 

As the size of the spatial directions is not bounded in any way, it has become 
customary to study spatial rather than temporal correlation functions. The spatial 
correlation functions are the static correlations in the equilibrium system, and hence 
are also of direct physical relevance. 

For the interpretation of the correlation functions measured in these simulations, 
one must discuss the spectral density functions underlying the correlations. While 
building models of the spectral densities, it is useful to keep certain constraints in 
mind. At finite temperatures the heat-bath provides a preferred frame of reference. 
As a result, the momentum representation of spectral functions is given in terms 
of po and \p\ separately, where p = (pi,P2,P3)- Dynamical modes are defined by 
poles in the complex po plane; and the movement of such poles with \p\ are the 
dispersion relations. In examining the spectral representation of the thermal and 
spatial direction correlators it becomes clear that different aspects of the spectral 
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function are important for each. However, information on the poles of spectral 
functions can be extracted from either of these correlators. 

The residue at a pole can also be used to extract effective couplings. These are 
usually defined as the integrals of discontinuities over cuts in the complex plane of 
either po or \p\. Note that on any lattice the spectral function does not have any 
cuts, but only a set of poles. However, effective couplings can still be extracted via 
sums over the residues at the poles. One can expects finite lattice spacing effects to 
come from those poles which develop into a cut in the continuum limit. 

2.2 Point Sources 

In the staggered discretization, local operators for currents carrying mesonic quan- 
tum numbers are written in terms of fermion fields, x( x )i as 

H(x) = (j) H (x)x(xo, x)x{xo, x). (2.3) 

The phase factors <pH project onto definite quantum numbers [[J, and are listed in 
Table [l], for the channels analyzed in this study. 

Since the meson operators, H( X^j , £1X6 composed of fermion operators, it is clear 
that even for a fixed external meson momentum, P or P, quarks with a spectrum 
of internal momenta will contribute to the correlation functions Gf and Gj* . The 
structure of the hadronic correlation functions is, therefore, more complicated than 
that of the quark correlator [0. 

This is clear when the correlation functions at O(g ) are explicitly written out. 
For the spatial correlation function the result is well known: 

24 ?V 3 

G?(x 3 ,P)/T 3 = —^Y.f H G q s (x 3 ,p)Gl(x 3 ,p>). (2.4) 



Table 1: Phases for the meson operators used for Gf . For Gf, x% should be re- 
placed by xq. Also listed are kinematic factors appearing in the O(g ) perturbative 
calculations of the meson correlation functions in eqs. ( |2.4| ) and fl2.6| ) . 



Channel 


(f) H (x) 


f H 

J even 


f H 
J odd 


S 


1 


m 2 — ujuj' 


1 


PS 




m 2 + UJUJ' 


1 


PV 


(-1) X1 + (-1) X2 + (-1)* 3 


3m 2 — Lou)' 


3 


V 




3m 2 + ujuj' 


3 
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The factors f H for even and odd sites are given in Table |T]. The quark propagator 
G q is given by: 



( smh[E s (x 3 - NJ2)] 
sinh(£ s iV CT /2)cosh£ s 

cosh[E s (x3-Ar a /2)] 



{ sinh(S s A^ CT /2)sinh2^ s 
The temporal correlation function is very similar: 



(x 3 odd) 
(x3 even) 



G?(x ,P)/T 3 



24iV? 



^f H G q t (x 0j p)GKx ,p') 



P 



with the quark propagator G\ given by 



G q t(x ,p) 



( cosh[E t (x - N T /2)} 
cosh(E t N T /2)coshE t 

swh[E t {x Q - N T /2)} 
[ cosh(E t N T /2)smh2E t 



(xq odd) 
(xq even) 



(2.5) 



(2.6) 



(2.7) 



The sums in eqns. (|2.4j ) and (|2.6| ) are over internal quark momenta. These run 
over po = (2n + l)irT and pi = 2nnT, with % = 1, 2, 3 and n = 0, ±1, • • • . The quark 
and antiquark momenta, p and p 1 respectively, sum to give the meson momentum 
P, i.e. = p^ + p' . We have used the abbreviation 



E 



t.s 



arcsinh(u; fiS 



{21 



along with the definitions 



u t = m l + sin2 (Pk) and 

k=l 



k=0 



(2.9) 



The variables E' and oj' are obtained by replacing p by p', and m q is the bare quark 
mass. 

Recall that the spectrum of po starts from the Matsubara frequency, Q = ttT, 
and reflects the antiperiodic boundary conditions usually imposed in the thermal 
direction on the fermion fields when inverting the Dirac equation. All other mo- 
mentum components start from zero. As a result, the screening mass at 0(g°) is 
independent of the meson operator H and given by H] 



Hh = %Ef n = 2arcsinh (mj + sin 2 fi) . 



(2.10) 
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Figure 1: Thermal direction correlation functions at O(g ) on lattices with N a — > oo 
and various N T values. The quark mass is m q = 0.01, and we define tT = xq/N t . 
The result for N T — > oo is obtained from eq. (|2.11|) . 

The decay of Gf at x 3 ^> 1/T is controlled by At shorter distances the 
higher modes also contribute and lead to a more rapid decrease of the correlation 
function. At very short distances, x% < 1/(2T), all terms in the sums in eq. ( |2.4| ) 
turn out to be important. 

The thermal direction correlation function, , does not yield extra information 
in perturbation theory. Since x® < 1/(2T), all terms in the sum in eq. (|2.6|) 
turn out to be important, and the correlation function falls very rapidly due to the 
contribution of terms with large E t and E' t . This is shown in Figure [l]. Also shown in 
the same figure is the result obtained letting first N a and then iV T — >• oo at fixed T. 



The formula? in eqs. (|2.6| ) and (|2.7| ) then reduce to the continuum expressions for the 
correlators constructed from non-interacting quarks. For zero external momentum, 
P = 0, and vanishing quark masses, these take on the simple form 

G?{t,P = Q) = — d y^-^co^\y{2Tr-l)) (2.11) 
ix Jo cosh y 

where y, t and T are related to the lattice variables via r = xqcl, tT = xq/N t and 
y = E t j (2T). Notice that, at short distances, one obtains G^(t, P = 0) ~ r~ 3 from 
eq. (|2.11|) . (Similar conclusions have been reached in the continuum, [fTOfl .) 

Although the temporal direction is too short to single out the lowest excitation 
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Figure 2: Scaled local masses at O(g ), m(r)r, versus the scaled separation tT = 
x /N T , obtained from eqs. ( |2.6| ) and (|2.7|) for N a — > oo for N T = 16 (dashed line) 
and 32 (dotted line). The full line shows the result for N T = oo (eq. ( |2.11|) ). The 
filled square shows the result for m(r = 0.25/T) on an 8 x 16 3 lattice. 

for T > 0, it is still instructive to study the behaviour of local masses, m(xo), 
obtained from G^ 1 . These are obtained by comparing the correlation functions on 
successive even (or odd) lattice sites, x ± 1. One then solves 

G?(x -1,0) cosh[m(xo)(x - 1 - N T /2)] 

Gf(x + 1, 0) cosh [m(x )(x + 1 - N T /2)] ' 1 ' ' 

for m(xo). The left hand side is obtained from either a measured or perturbatively 
calculated correlation function. One then assumes on the right hand side that the 
correlation function can be described by an ordinary meson correlator consisting of 
a single cosh. (For the simple case of non-interacting quarks we know, of course, 
that we should rather use the sum over squares of quark correlators applicable to 
the given source and sink.) A feeling for this local mass may be obtained by noting 
that in the limit iV r — > oo, at fixed T (i.e., a — *■ 0), eq. ( p.ll ) can be used to show 
that 

m(x )x = m(r)r ~ 3 (iV r > x > 0). (2.13) 

Analytic results for the local masses at O(g ) in perturbation theory are shown 
in Figure [| Note that, in practice, m(r)r remains close to 3 for all values of tT. 
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The rise of the local meson mass at short distances, when extracted in this fashion, 
reflects the increasing importance of higher momentum excitations. We note that 
even for the largest possible separation x = N T /2, or equivalently r = 1/(2T), 
the local mass remains large in units of the temperature, m/T « 6. The result for 
t = 1/ (4T), obtained on an 8 x 16 3 lattice, is also shown in the figure. This is the 
only value of r which we will be able to study numerically. 

As a concluding remark to this section, let us point out a distinctive feature of a 
meson propagator consisting of non-interacting quarks, GJ*(xq), namely the strong 
oscillation between even and odd sites visible in Figure |l|. This is seen for any 
non-zero lattice spacing (finite N T ), and persists in higher orders of perturbation 
theory. Such oscillatory behaviour is also seen in the spatial correlation functions 
in perturbation theory. In the low-temperature phase of QCD, however, where one 
has a genuine bound state, the correlator does not show such a behaviour. 

2.3 Wall Sources 

From the previous discussion we conclude that the analysis of the correlation func- 
tions (and especially the temporal correlation functions) between local hadronic 
operators is complicated, owing to the momentum sums in eqs. ( p.4|) and ( |2.6|) . A 
determination of the low momentum excitations through local masses or fits with 
a few exponentials seems to be difficult, as all higher energy levels, coming from 
higher quark momenta, contribute. 

The contributions from higher quark momenta can be suppressed by the judicious 
choice of a wall source, which creates quarks with only the lowest allowed momentum. 
For correlations measured in the thermal direction, the boundary conditions on the 
transverse slice are periodic. In this case one can project onto zero quark momentum 
via the meson sources 

g 

H wall (x , x) = — J2 0h(x)x(xo, x + e)x(x , x + e). (2.14) 
° e,e' 

even 

For correlators measured in the spatial direction, the antiperiodic boundary con- 
dition in the thermal direction means that one can not project onto zero momentum. 
Instead a projection onto the Matsubara momentum f2 = (ttT, 0, 0) must be used: 

H Q (x 3 ) = <Ph{x)x{x3, x + e)x{x 3 , x + e!) exp (i(e - e!) ■ o) . (2.15) 

even 

The correlation functions are measured between a wall source and a point sink, with 
a sum over all sinks in the plane transverse to the direction of propagation. 
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The 0(g°) perturbative calculation in eqs. fl2.4Q and (|2.6| ) is simplified if the only 



contribution is from quarks with momentum p = or p = ft. For example, the 
temporal correlator becomes the square of a single quark propagator, 

Gf' wall (x ,0)/T 3 = 3(N 3 T /N 3 a )f H G!(x ,0)G!(x ,0). (2.16) 

The biggest gain in using extended sources is obtained for the thermal direction 
correlation functions. At O(g ) in perturbation theory the masses are then given by 



m 



wal1 = 2arcsinh (m q ) . (2.17) 



Local masses can then be extracted by using eq. (|2.12|) . We note that at O(g ) the 



meson correlator now has the form cosh 2 (m q (xo — N T /2)), which is well fitted by a 
cosh(2m q (xo — N T /2)) over most of the temporal range. 

At O(g ) in perturbation theory the local masses thus measure the mass of the 
lowest excitation. At higher orders this is modified due to diagrams which can be 
separated into thermal corrections to each of the quark propagators and interactions 
between the propagating quark and antiquark. Each of these effects begins at O(gT). 
If the effective interaction between the two quarks can be neglected, as we expect for 
the V channel, then the effect of interactions can be wholly subsumed into replacing 
the quark mass m q in eq. (|2.17|) by the thermal quark mass, m^. This was extracted 
from the temporal quark correlator in Q using the same configurations used here. 
It is not unrealistic to hope that such a result may hold beyond perturbation theory 
as well. 



2.4 Varying Boundary Conditions 

In our attempt to clarify the nature of the excitations in the plasma we have also 
varied the temporal boundary condition on the valence fermion fields from which 
the meson operators are constructed. The boundary conditions of the sea quarks, 
driving the generation of the gauge field ensemble, have been left untouched. This 
should give a direct view of the nature of excitations in the plasma. 

If the lowest excitation in a channel of fixed quantum numbers consists of a me- 
son, then the spectral representation contains a pole. Thus, we are assured that this 
state will be seen in the correlation of generic operators with these quantum num- 
bers. In particular, whether we use periodic or anti-periodic boundary conditions on 
the fermion source, the lowest mass obtained out of the correlation functions must 
be the same. The correlation functions should thus show little dependence on the 
choice of boundary conditions. 
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However, if a given quantum number can only be obtained by exchanging more 
than one particle, then the spectral representation contains a cut. If the exchanged 
particles are fermions, then the contribution to the correlation function depends on 
the boundary conditions. In fact, at O(g ), it is easy to see that eqs. ( |2.4| ) and Q2.5|), 
which define the spatial correlation functions, depend on the temporal boundary 
conditions via the spectrum of po alone. For periodic boundary conditions in the 
thermal direction, the allowed values are p Q = 2nirT with n = {0, 1, • • •}, which 
yields a screening mass equal to 

H = 2Ef n = 2arcsinh(m q ). (2.18) 

The qualitative change in going from anti-periodic to periodic boundary conditions 
in the thermal direction is thus expected to be significant, and should show up in 
the correlation functions. Interactions will modify this in the same way as discussed 
above. 



2.5 Effective inter-quark Couplings 

One may use the meson propagators to calculate an effective four fermion coupling 
as used in, for example, the Nambu-Jona-Lasinio model |]llj. The propagator at 
zero four- momentum is equal to the generalised susceptibility 

X,,p = G v , p (0). (2.19) 

Letting G° n (0) = x1, P be the susceptibility for non- interacting quarks in the 
pion and rho channels respectively, one obtains from the Dyson equation the result: 

X,, P = ,_ X }'"o » ( 2 - 2 °) 

where g WlP /2 is the four fermion coupling in the effective Lagrangian, following the 
convention of ||11|| . 

Solving for g n<p one obtains: 

9n, P T 2 = ( - Q -) , (2.21) 

N r \Xn,p X^p) 

where the susceptibility extrapolated to zero quark mass is used. 

°The corresponding formula, eq. 19, given in 0], contains a minor error in the normalization. 
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3 Results 



In this section we describe the results of our measurements on configurations gener- 
ated by the MT C collaboration on lattices of size 8 x 16 3 with 4 flavours of dynamical 
fermions having a bare mass m q a = 0.01. Recall that the phase transition was ob- 
served at a coupling of /3 — 5.15(5) fll2| . At this coupling we have two sets of 
configurations. The set labelled B corresponds to a run which started in the low- 
temperature (chiral symmetry broken) phase and remained there. The set labelled S 
was started in the high-temperature (chirally symmetric) phase and did not tunnel 
into the other phase. The couplings (3 = 5.1 and 5.2 cover the temperature range 
T/T c = 1 ± 0.2. 

Results are presented for both point sources and wall sources. For the wall 
sources we have fixed the configurations to the Coulomb gauge in the hyperplane 
containing the source, transverse to the direction of propagation. 

We have around twenty configurations at each value of the coupling, with four 
sources per configuration for point sources, and one for wall sources. In calculating 
the errors for the point sources we first blocked the four sources. The errors quoted 
reflect the statistical fluctuations alone. 

3.1 Thermal direction point source correlators 

The measured values of Gf are shown in Figure |3|. Note that the pseudoscalar 
correlator for T < T c is very well described by a hyperbolic cosine. A changeover to 
the oscillatory behaviour characteristic of free fermions, as discussed in the previous 
section, is visible in our thermal correlation functions as the critical coupling, (3 C = 
5.15, is crossed. This change in behaviour is most noticeable when wall sources are 
used. Note also that the oscillation between even and odd sites is enhanced in the 
vector channel. This is in accordance with the perturbative calculation, and is due 
to the fact that f v > f PS (see Table 0). 

In Figure [|we show the local masses at distance r = 0.25/T extracted from the 
V and PS correlation functions using eq. (|2.12|) , (and hence the implied ansatz,) in 
the thermal direction. These are also listed in Table |2|. We find that the local mass 
in the V channel is close to the perturbative value, my = 9.4/T (see Figure 0). 
However, the local mass in the PS channel approaches this value rather slowly with 
increasing temperature. Such a behaviour is very similar to that of the screening 
masses extracted from spatial correlation functions. 

These local masses are similar in magnitude to the screening masses. This is 
accidental. It is due to the fact that for N T = 8, N a = 16 we can only extract a local 



11 




Figure 3: The thermal direction correlators for PS (filled circles) and V (filled 
squares) at (3 = 5.1 (a), 5.3 (b) and 6.5 (c), shown as a function of tT. In (c) 
the lowest order results for PS (solid line) and V (dashed line) correlators are also 
shown. 

temporal mass at tT =1/4. In a free fermion theory on this size lattice we find 
that the local mass at this distance is about 10T. This just happens to be close to 
the meson screening mass in free fermion theory on this size lattice, 9T. When the 
lattice size is changed, this accidental concordance is removed (see eq. ( |2.13| )). 



Table 2: Local masses from thermal direction V and PS correlation functions at 
distance r = 1/(4T). Results obtained from point and wall source operators are 
shown. For comparison we also list twice the effective quark mass J7j. 



p 


point 
m PS 


point 

m v 


/Krl wall 
m PS 


, rrl wall 

Illy 


2mf 


5.1 


0.24(5) 


1.5(2) 






0.9(1) 


5.15(B) 


0.36(7) 


1.2(1) 






0.6(1) 


5.15(S) 


0.46(6) 


1.19(9) 


0.29(1) 


0.56(1) 


0.42(6) 


5.2 


0.51(6) 


1.08(8) 


0.31(1) 


0.493(5) 


0.38(4) 


5.3 


0.64(4) 


1.13(6) 


0.303(6) 


0.401(5) 


0.19(9) 


6.5 


0.98(2) 


1.13(3) 


0.0838(4) 


0.106(2) 


0.044(8) 
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Figure 4: Local masses m H ° at distance tT = 1/4 from thermal direction 
V (squares) and PS (circles) correlation functions using point sources, shown as 
a function of (3. Also shown is the corresponding result computed in a theory of 
non- interacting fermions on an 8 x 16 3 lattice at the same r (full line). 
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Figure 5: Local masses m H from correlators constructed with wall sources 

in the PS (circles) and V (squares) channels in the high temperature phase. Also 
shown is 2m^ s (triangles) from 0. Lines have been drawn to guide the eye. 
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3.2 Thermal direction wall source correlators 



Local masses have similarly been extracted from the temporal wall source correla- 
tors. The results obtained for m wa11 are also collected in Table 0. We note that 
these masses are indeed much smaller than the local masses obtained using point 
sources. The projection onto the lowest momentum excitation with these correlation 
functions thus seems to be rather efficient. 

The relation in eq. fl2.17|) seems to hold, at least qualitatively, for both the vector 



and the pseudo-scalar channels, when the effective quark mass given in |7] is 
used for the quark mass. The masses in the PS channel close to (3 C are smaller 
than those in the vector channel, indicating strong residual interactions between 
fermions. The PS masses then approach the mass in the vector channel as the 
coupling increases. 

We also note that the vector mass remains larger than 2m^, indicating the 
importance of the remaining inter-quark interactions. One may try to parametrize 
these residual interactions in a potential model relating the difference in the pion and 
rho masses to different spin-spin interactions in these quantum number channels. 
The masses then receive contributions from the effective quark masses, as well as the 
scalar (-E sca iar) and spin dependent (E spin ) part of the quark-anti-quark potential: 

3 

ITLps = 2mf + -E SC alar - ^spin (3.1) 

1 

m v = 2mf + E Bcalax + -£ spin . (3.2) 

We find that the meson masses can then be parametrized by a scalar term that is 
consistent with zero, and a spin term approximately equal to the effective quark 
mass, and with the same temperature dependence. 

3.3 Anti-Periodic and Periodic Boundary Conditions 

In Figure |6| (a) and (b) Gf s and are shown for T < T c for both periodic (per) 
(a) and anti-periodic (aper) (b) boundary conditions in the thermal direction. The 
correlation functions are unaffected by this change, and the screening masses, shown 
in Table |3], do not change within errors. Thus, these masses reflect bosonic poles in 
the spectral function. 

In contrast, there is a remarkable difference between the correlation functions 
obtained for periodic and anti-periodic boundary conditions for T > T c , as shown 
in Figure |6| (c) and (d). The correlation functions become much flatter when the 
boundary conditions are periodic, and the screening mass drops in both the PS and 
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Figure 6: in the PS (circles) and V (squares) channels. Figures (a) and (b) 
correspond to (3 = 5.1, while figures (c) and (d) correspond to (3 — 5.3. The valence 
quarks have anti-periodic boundary conditions in the thermal direction for (a) and 
(c), with (b) and (d) having periodic boundary conditions in the thermal direction. 
The masses shown were obtained from fits. 



V channels. Varying the boundary conditions thus suggests the absence of genuine 
mesons in the high temperature phase of QCD. 

It can be seen that the vector screening mass /Zy er closely satisfies the relation 
in eq. ( 2.18 ), when the quark mass is taken to be the non-perturbative quark mass 
measured at the same coupling [|7|. Thus, the effects of interactions, in this an- 
gular momentum channel, can be almost entirely lumped into the effective quark 
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mass. (This relation also holds below the transition, and gives a definition of the 
constituent quark mass.) As is already known ||, this simplification does not hold 
in the PS channel, and an effective interaction between quarks remains. The value 
of /Xpg is accordingly somewhat smaller than /Zy er . 

In Table |3| and Figure [F| results for screening masses from the wall source op- 
erators (eq. fl2.15|) ) are presented. These are a little smaller than those obtained 



from point source operators, which indicates that the spatial direction is not large 
enough to eliminate higher terms in the sum in eq. (|2.4j ). Presumably, slightly larger 
spatial lattices would be required for this. Experience from j|] shows that spatial 
sizes N a « 4N T generally suffice to eliminate the effect of the higher modes. 

The values of /J^ 1 are similar to the values of local masses extracted from wall 
source operators in the thermal direction. This is indeed to be expected in pertur- 
bation theory. At O(g ) these two masses should be the same, while differences can 
arise at higher orders in g. Thus in the limit m q — > 0, both these quantities are 
sensitive to O(gT) thermal corrections. The screening masses in Figure |7| may be 
compared with the wall source results shown in Figure [5[ 

3.4 Effective inter-quark Couplings 

The effective four fermion coupling, defined in eq. ( j2.21| ), is presented in table |j. 
Notice that the coupling in the pion channel is about four times stronger than that 



Table 3: Screening masses in the PS and V channels for both anti-periodic (aper) 
and periodic (per) boundary conditions in the thermal direction. In the chirally 
symmetric phase results obtained with wall source operators (W) are also quoted. 
The screening mass expected for this size lattice, using non-interacting quarks, is 
2ttT = 1.2 for point sources and 2irT = 1.0 for wall sources. 
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Figure 7: Screening masses extracted from spatial correlation functions using point 
sources with periodic temporal boundary conditions. Shown is the rho (squares) 
and pion (circles) screening mass, extracted in both cases from fits to the correlator 
on even sites. Also shown is twice the effective quark mass (triangles) determined 
on the same gauge field configurations in the Landau gauge 0. Lines are drawn to 
guide the eye. 



in the rho channel, supporting the hypothesis that the difference between the pion 
and rho screening masses lies in different interaction strengths between the quarks 
in the two channels. The numbers obtained at the transition are also comparable 
with those obtained in the quenched approximation Q. 

Taking the transition temperature to be T c = 0.14GeV [O the following phys- 



Table 4: Effective four fermion couplings for the pion (g n ) and rho (g p ) in the limit of 
zero quark mass. The susceptibilities for non-interacting quarks are: l/x„ = 0.5478 
and l/x° p = 0.3680. 
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0.005(8) 
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ical values for the couplings in the chirally symmetric phase at the transition are 
obtained: g n = 0.41GeV~ 2 and g p = 0.094GeV -2 . These numbers may be compared 
with the values quoted in obtained from fits to the experimental meson data 
at T = 0: g n = 4.90GeV~ 2 and g p = 3.25GeV~ 2 . The couplings at high tempera- 
ture are well below the critical coupling at which the Nambu-Jona-Lasinio model 
first shows chiral symmetry breaking, and provide a further indication that neither 
channel has a low lying bound state. 



4 Conclusions 

The pion and rho correlators are, above the phase transition, sensitive to both the 
boundary conditions and the type of source used. This is not seen below the phase 
transition. Since the changes in the correlator are of the form one expects if unbound 
fermions play a direct role in the spectral function, this provides evidence for the 
existence of a two fermion cut dominating the spectral function. 

Above the phase transition the screening masses in both the PS and vector 
channels are consistent with twice the effective quark mass plus some residual inter- 
actions. As a measure of this interaction, and hence as a summary of the relevant 
physics of the system, we extracted an effective four fermion coupling. This was 
four times stronger for the PS channel than it was for the vector channel, and an 
order of magnitude smaller than the couplings used in Nambu-Jona-Lasinio models 
at zero temperature. 

The structure of correlators in the vector channel above T c generally agrees quite 
well with the behaviour expected from leading order perturbation theory; however, 
this is not the case for the pseudo- vector channel below 2T C . Here the correlators 
and masses are seen to approach the perturbation limit rather slowly. 

One is left with a consistent picture of a plasma phase consisting of deconfined, 
but strongly interacting quarks and gluons in the temperature range from T c to 2T C . 
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